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Popular Summary 

TRMM (Tropical Rainfall Measurement Mission) offers a unique opportunity to 
improve the understanding of tropical meteorology and to evaluate the impact of rainfall 
data on tropical weather forecasts. Early studies demonstrated that assimilation of the 
TRMM microwave imager (TMI) derived rainfall data into large-scale global model is 
beneficial for the analysis of atmospheric general circulation (Hou et al. 2000) and also 
consequently can have significant impact on mesoscale forecast of Supertyphoon Paka 
(1997) (Pu et al. 2002). Instead of direct assimilation of rainfall data into the global 
model as in Hou et al. (2000), this paper evaluates the impact of TMI rainfall on 
mesoscale forecast via the direct assimilation of the TMI-derived rainfall rate into the 
mesoscale regional model itself using a four-dimensional variational data assimilation 
(4DVAR) technique. 

Sensitivity studies are performed on the assimilation of TMI-derived rainfall data 
into a mesoscale model using a four-dimensional variational data assimilation (4DVAR) 
technique. A series of numerical experiments is conducted to evaluate the impact of TMI 
rainfall data on the numerical simulation of Hurricane Bonnie (1998). The results indicate 
that rainfall data assimilation is sensitive to the error characteristics of the data and the 
inclusion of physics in the adjoint and forward models. In addition, assimilating the 
rainfall data alone is helpful for producing a more realistic eye and rain bands in the 
hurricane but does not ensure improvements in hurricane intensity forecasts. Further 
study indicated that it is necessary to incorporate TMI rainfall data together with other 
types of data such as wind data into the model, in which case the inclusion of the rainfall 
data further improves the intensity forecast of the hurricane. This implies that proper 
constraints may be needed for rainfall assimilation. 

In addition, as the coverage of TMI data is very limited for the regional 
applications, the use of merged multi-satellite data or future Global Precipitation Mission 
(GPM) data products will be necessary for future applications. 
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ABSTRACT 


Sensitivity studies are performed on the assimilation of TRMM (Tropical 
Rainfall Measurement Mission) Microwave Ima ger (TMI) derived rainfall data into a 
mesoscale model using a four-dimensional variational data assimilation (4DVAR) 
technique. A series of numerical experiments is conducted to evaluate the impact of TMI 
rainfall data on the numerical simulation of Hurricane Bonnie (1998). The results indicate 
that rainfall data assimilation is sensitive to the error characteristics of the data and the 
inclusion of physics in the adjoint and forward models. In addition, ass imilatin g the 
rainfall data alone is helpful for producing a more realistic eye and rain bands in the 
hurricane but does not ensure improvements in hurricane intensity forecasts. Further 
study indicated that it is necessary to incorporate TMI rainfall data together with other 
types of data such as wind data into the model, in which case the inclusion of the rainfall 
data further improves the intensity forecast of the hurricane. This implies that proper 
constraints may be needed for rainfall assimilation. 
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1. Introduction 


The Tropical Rainfall Measuring Mission CTRMM) is a joint Japan-U.S. project 
to measure rainfall over the global tropics. With the world’s first precipitation radar, the 
TRMM satellite has provided the first detailed and comprehensive dataset on the four- 
dimensional distribution of rainfall and latent heating over the tropics (between 35° N and 
35° S). TRMM offers a unique opportunity to improve the understanding of tropical 
meteorology and to evaluate the impact of rainfall data on tropical weather forecasts. 
Early studies have demonstrated that assimilation of TRMM microwave imager (TMI) 
derived rainfall data into large-scale global models is beneficial for the analysis of the 
atmospheric general circulation (Hou et al. 2000) and also consequently can have 
significant impact on mesoscale forecasts [e.g. Supertyphoon Paka ini 997 (Pu et al. 
2002]. Instead of direct assimilation of rainfall data into the global model as in Hou et al. 
(2000) and Pu et al. (2002), this paper evaluates the impact of TMI rainfall on mesoscale 
forecasts via the direct assimilation of TMI-derived rainfall rates into the mesoscale 
regional model itself using a four-dimensional variational data assimilation (4DVAR) 
technique. 

The description of the model, the 4DVAR system, the hurricane case, and the TMI 
data is addressed in section 2 and 3, respectively. Detailed numerical results, including 
sensitivity studies and forecast impacts, are given in section 4. A summary and discussion 
are given in section 5. 

2. Description of the mesoscale model and 4DVAR system 


2 


The Penn State University/National Center for Atmospheric Research 
(PSU/NCAR) mesoscale forecast model (MM5) and its adjoint system are used in this 
study. The MM5 is a limited-area, non-hydrostatic primitive equation model with 
multiple options for various physical parameterization schemes (Dudhia 1993; Grell et al. 
1995). The model employs a terrain-following a vertical coordinate, where cs is defined 
as o- (p - PtopV (p<& “Ptop)’ P pressure, and p s f C and p top are the pressures at the 

surface and model top, respectively. Physics options used for the forecast model in this 
study include the Grell cumulus parameterization, a simple ice microphysics scheme 
(Dudhia 1993), the Blackadar high-resolution planetary boundary layer parameterization 
scheme (Blackadar, 1976, 1979; Zhang and Anthes 1982), and a cloud atmospheric 
radiation scheme (Dudhia 1993). The land surface temperature is predicted using surface 
energy budget equations as described in Grell et al. (1995). For a more detailed 
description of MM5, see Dudhia (1993) and Grell et al. (1995). 

The MM5 adjoint modeling system (Zou et al. 1998) is employed in the data 
assimilation experiments. For the variational data assimilation system, the physics options 
in the adjoint model are the Grell cumulus parameterization, a simple ice microphysics 
scheme (Dudhia 1989), and the Blackadar high-resolution planetary boundary layer 
parameterization scheme. Application of the MM5 adjoint model to a variety of mesoscale 
weather systems has been demonstrated in papers by Kuo et al. (1996) and Zou and Kuo 
(1996). 

In general, a 4DVAR system tries to minimize the following cost function: 

J(x.)='ZJ k +J l , ( 1 ) 
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where x is the analysis variable and the subscript “0” denotes the initial state. J b is the 
background term, and is the contribution to the cost function from an individual type of 
data. The subscript k denotes the type of data and m is the total number of available data 
types. For example, the contributions from one arbitrary type of observation can be 
described as follows: 

.(M.iX'V-oJW . ( 2 ) 

1=0^ 

where O is the observation data, “i” denotes the “ir/z” time step for integration of the non- 
linear forecast model M, at which the observations are available, and i € (0, A) > while A 

is the length of the assimilation window. W is a weighting factor that depends on the 
statistical error characteristics of the observational data. H is a so-called observation 
operator (possibly no-linear), which transfers the grid-space model variable x to the 
observational type. In order to minimize the cost function, the adjoint of the tangent 
linear model of the nonlinear forecast model is required (Talagrand 1987). 

The effectiveness of the 4DVAR technique for assimilation of precipitation 
observations has been addressed by Zou and Kuo (1996) and Zupanski and Messinger 
(1995). However, this paper will further investigate strategies on the assimilation of TMI 
data particularly. 

3. Hurricane Case and TMI rainfall data 

The TRMM Microwave Imager (TMI) is one of several TRMM satellite sensors. 
The TMI measures the horizontal distribution of rainfall by receiving microwaves 
emitted or scattered by raindrops and ice particles in five microwave channels. At NASA 
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Goddard Space Flight Center, the TMI microwave radiances are used to retrieve surface 
rainfall rate information via the Goddard Profiling (GPROF) algorithm. The basis for the 
rainfall retrieval algorithm is the Bayesian technique described in Kummerow et al. 
(1996) and Olson et al. (1996, 1999). In order to evaluate the impact of TMI rainfall data 
on mesoscale forecasts, these retrieved surface rainfall data are assimilated into the MM5 
model. 

The TMI footprints usually cover the global tropics (35° N-35° S) in a 24-hour 
period. However, during a typical mesoscale analysis period (usually a 3 h or 6 h period) 
and for a specific regional domain, there are only limited TMI observations available. In 
most cases, the TMI samples only about twice a day for a certain region, and the time 
interval between the swathes may exceed 6 h. Considering the data availability during a 
6h analysis cycle. Hurricane Bonnie (1998) was selected from several storm cases to 
perform the sensitivity studies in this paper. 

There were two TRMM swathes that passed over Bonnie in the Atlantic Ocean with 
a time interval of about 6 h. The two overpasses were around 1139 UTC 22 August 1998 
and 1807 UTC 22 August 1998 (Fig.l), respectively. At the time, Bonnie was a category 1 
hurricane based on the Saffir-Simpson intensity scale, having recently developed from a 
tropical storm. 1200 UTC 22 August 1998 was selected as the initial time for the 
experiments. A 6 h data assimilation window was set for the period 1200 UTC - 1800 UTC 
August 1998. 

Based on TMI-derived surface rainfall data being typically defined as an hourly 
“rain rate” and the actual data availability, the TMI-derived rainfall was treated as “hourly 


( 3 ) 
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averaged rainfall”. Therefore, Eq.(2) can be written as follows: 




f 1 17 


= 0 ,A 


^2 CRv-RR , 

1 j-0 


f 1 .5 




^icRi-RR, 

1 j=0 


where T is the averaging time period (i.e., one hour), RR the retrieved TMI rain rate, and 
CR the model-generated rainfall in the time step. 


4. Numerical Experiments and Results 

The data assimilation experiments were conducted at 36 km resolution. The 
model domain is show in Fig.l. 


4.1 Sensitivity studies 

Numerical experiments were conducted to test strategies for assimilating TMI- 
retrieved rainfall rates. Two groups of sensitivity studies were performed to test the 
sensitivity of specification of the error characteristics and the inclusion of physics in the 
adjoint model to the TMI rainfall data assimilation. Table 1 lists the experimental 
configuration for all numerical experiments in this paper. 

a. Sensitivity of the specified rainfall data error characteristics 

The W factor in Eq.2 is a weighting factor that depends on the statistical error 
characteristics of the observations. To some extent, this factor represents how much the 
4DVAR system would “trust” the observations. Because the correlations between the 
observations are usually unknown or difficult to define, the W matrix is often defined as a 
diagonal matrix. Three experiments were conducted with the specification of the error 
characteristics as follows: in Experiment 1, W was set up as a unified number and 
defined as the inversion of variances based on all available data (i.e., unified W); in 
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Experiment 2, a 20% error was assumed for the retrieved rainfall rate, and W was defined 
as the inversion of the error variances (i.e., 20% error); and in Experiment 3, the error 
characteristics were specified following Bauer et al. (2002) and Olson (personal 
communication) as: 


over ocean 

CTs =1.357 R 0 ' 7 

over land 

CTs =2.516 R 0 - 558 


where R is the retrieved rain rate and ct s the standard deviation of the rain rate. The a s 
was obtained from a large sample of retrieved rainfall rate datasets. 

Fig.2 shows the variation of the cost-function with the number of iterations. The 
definition of W mainly impacts the convergence of the minimization in terms of both 
cost-function reduction and speed of convergence. As a consequence, W effects how 
much information can be gained from the observations. The results show that it is 
obviously advantageous to use the error specification suggested by Bauer et al. (2002 ) 
and Olson (personal comm.). 

b. Inclusion of physics in the adjoint 

For a common forecast model (forward model), the cumulus parameterization and 
microphysical processes usually help the model to produce a better rainfall forecast. 
However, due to the difficulties in deriving an adjoint model for the physics package, in 
some previous studies (e.g., Zou and Kuo 1996), not all of the physics processes were 
included in the adjoint model. In order to test the impact of including physics in both the 
forward and adjoint models on data assimilation results, the following three experiments 
were conducted: in Experiment 4, both forward and adjoint models include cumulus 
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parameterization but not microphysics; in Experiment 5, both cumulus parameterization 
and microphysical processes are included in the forward and adjoint models; and in 
Experiment 6, neither cumulus parameterization nor microphysical process are included 
in adjoint and forward models. For all Experiment 4-6, the specified error characteristics 
followed those in Experiment 3. 

Fig.3 shows the variation of the cost-function with the number of iterations. As 
expected, including all of the physics packages has the largest benefit in terms of both 
convergence and assimilation results. 

4.2. Impact on forecasts 

Figure 4 shows rainfall rates at the end of the data assimilation (6-h forecast) from 
Experiment 2 compared with the control experiment (CTRL) where rainfall data were not 
assimilated. Obviously, assimilating the rainfall data helps the model to produce a more 
realistic eye and rain bands in the hurricane. The results are quite encouraging. Twenty- 
four hour forecasts were conducted to test the impact of the rainfall data on the hurricane 
intensity forecast. Unfortunately, there was no improvement in the consequent forecasts. 
The forecasted track and intensity (in terms of maximum surface wind and minimum sea 
level pressure) are almost the same in cases both with and without the TMI rainfall data 
assimilation (figure not shown), indicating the impact of rainfall data on consequent 
forecasts is almost negligible. This is not consistent with previous results (Pu et al. 2002). 

In order to explore additional strategies for rainfall data assimilation, an additional 
experiment was performed to assimilate rainfall data along with other data sets. As other 
conventional data is unavailable for this case, bogus vortex wind information is 
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introduced into the assimilation process. Following Pu and Braun (2001), the bogus wind 
data was derived from the gradient balance equations as follows: 

C" (r) = [AB(p. - Pc )exp(- A/r B )/ pr* f (4) 

where \4 >osus is the gradient surface wind at radius r, p the air density (assumed constant 

at 1 . 15 g m' 3 ), p c the central pressure and p n the ambient pressure (theoretically at infinite 
radius, however, here taken from representative values in the hurricane environment). 
The scaling parameters A and B are defined by maximum wind information. By setting 
dVj dr = 0 , the radius of maximum surface wind (RMW) is R m = A ],B , and substitution 

back into (4) gives the maximum wind speed, V m = C(p n - p c ) n , where C = (B/ pe) n and 
e is the base of the natural logarithm. 

Based upon the best available estimates (according to the report from the 
Hurricane Research Division/AOML/NOAA), the parameters defining the bogus vortex 
are given by p c =980hPa centered at (22.3°N, 69.8°W), p n =1012hPa, V m =38.6 m s' 1 , and 
R m =120 km. The bogus wind information extends out to a radius of 350km. The surface 
wind is then extended into the vertical with a vertical profile following Pu and Braun 
(2001). For the experiments, the specified wind information is assimilated every 10 
minutes wit hin the first 30 minutes. 

Two sets of numerical experiments were performed in a 6-h assimilation window. 
In the first experiment (Experiment 7), only bogus wind information was assimilated into 
the model. In the second experiment (Experiment 8), the rainfall data are incorporated 
along with the bogus wind information. 

Figures 5 and 6 show the rainfall rates at the end of the data assimilation (e.g., 6- 
h forecasts) for Experiment 7 and 8, respectively. Compared with the TMI observations 
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(Fig.lb), the rainfall patterns in Fig. 8 are closer to the observed rainfall structure, 
indicating that assimilation of rainfall data further improves the asymmetric hurricane 
rainfall structure. 

Further comparison is illustrated by histograms of the probability density function 
(PDF) of 1-h rainfall amounts at the end of the data assimilation (Fig.7). The figure 
shows that rainfall is generally underestimated in the case with rainfall data assimilation 
but overestimated in the case without rainfall data assimilation. However, the spectral 
distribution of rainfall rates is relatively narrow, with only one peak in the case with 
rainfall data assimilation. When the rainfall data are assimilated into the model, the 
spectrum of rainfall rates becomes broad with multiple peaks; the heavy rainfall rates 
(~1 O mm/ hr) are also better represented compared to the case without rainfall 
assimilation. 

Figure 8 shows the time variation of the forecast hurricane intensity in terms of 
the maximum surface wind and minimum sea-level pressure for the subsequent 24 h 
forecasts. It indicates that rainfall data assimilation is not only helpful for producing 
better vortex rainbands but also improves subsequent forecasts. The positive impact 
shown in this group of experiments indicates that it may be necessary to incorporate 
TRMM r ainf all data together with other types of data such as wind data in order to 
further improve the intensity forecasts for hurricanes. This implies that some constraints 
may be needed for rainfall assimilation. 

5. Summary and discussion 
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The main conclusions from numerical experiments with TMI rainfall assimilation 
in this study are: 

• Rainfall data assimilation is sensitive to the error characteristics of the data and 
the inclusion of physics in the forward and adjoint models, suggesting that it is 
necessary to use the full physics model in rainfall data assimilation and to take 
into account the error characteristic of the data; 

• Ass imilating the r ainfall data alone produces a more realistic eye and rain 
bands in the hurricane but does not ensure improvements in hurricane intensity 
forecasts. Numerical results indicate that it is necessary to incorporate TRMM 
rainfall data together with other types of data such as wind data into the model, 
in which case the inclusion of the r ainfall data will further improve the intensity 
forecast of the hurricane. This implies that some constraints may be needed for 
r ainfall as similati on. 

In addition to the TMI data, a similar experiment was performed assimilating 
surface rainfall data derived (from the same algorithm, i.e. GPROF) from TRMM 
precipitation radar (PR) for Hurricane Bonnie for the same assimilation window. 
Fortunately, the PR swathes overlap the TMI swathes in both time and space except that 
the PR swathes are much narrower than the TMI swathes (i.e., they cover the one third of 
the TMI swathes, figure not shown). A data assimilation experiment similar to 
Experiment 3 was performed with the PR data. The results are very similar to those from 
the TMI data assimilation, suggesting that the PR can also provide useful data sets for 
rainfall assimilation. 
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Future studies will be conducted to further confirm the above conclusions, and to 
explore the possibility of incorporating TMI and PR rainfall with other conventional and 
satellite data to improve mesoscale precipitation and storm forecasting. On the other 
hand, the coverage of both TMI and PR data is very limited for regional applications. 
Therefore, the use of merged multi-satellite data (e.g., Huffman et al. 2001; Huffman and 
Bolvin 2003) will be another alternative option for future study. However, since the 
Global Precipitation Mission (GPM) is under preparation, a large benefit to NWP could 
be obtained from that global precipitation data. 
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Table 1. Experimental Design 


Experiment 

Number 

Definition 
of W Term 
for rainfall 

Physics 

Included 

Bogus 

Vortex 

Rainfall 

Assimilation 

CTRL 

No 

Grell cumulus scheme 

No 

No 

1 

Unified 

Dudhia microphysics 
Grell cumulus scheme 

No 

Yes 

2 

20% error 

Dudhia microphysics 
Grell cumulus scheme 

No 

Yes 

3 

Bauer 

Dudhia microphysics 
Grell cumulus scheme 

No 

Yes 

4 

Bauer 

Dudhia microphysics 
Grell cumulus scheme 

No 

Yes 

5 

Bauer 

Grell cumulus scheme 

No 

Yes 

6 

Bauer 

Dudhia microphysics 
None 

No 

Yes 

7 

No 

Grell cumulus scheme 

Yes 

No 

8 

Bauer 

Dudhia microphysics 
Grell cumulus scheme 
Dudhia microphysics 

Yes 

Yes 
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Figure Caption 


Fig.l Rain rates (mm/hr) for two TMI swathes that passed over the Hurricane Bonnie 
(1998) around a) 1139 UTC 22 August 1998 and b) 1807 UTC 22 August 1998. 

Fig.2 Variation of normalized cost-function with iteration number for Experiment 1 
(dotted line). Experiment 2(dashed line) and Experiment 3 (solid line). 

Fig.3 Same as Fig.2, but for Experiment 4 (dotted line). Experiment 5 (solid line) and 
Experiment 6 (dashed line). 

Fig.4 Hourly accumulated rainfall rate (mm/h) at the end of data assimilation (6h) for a) the 
control run (CTRL), b) Experiment 5 and c) the differences between a) and b). 

Fig.5 Hourly accumulated r ainf all rate (mm/h) at the end of data assimilation for 
assimilation of bogus wind only (Experiment 7). 

Fig.6 Same as Fig.5, but both bogus wind and rainfall data assimilated (Experiment 8). 

Fig.7 Histograms of probability density functions of 1-h rainfall amounts at the end of data 
as similati on for Experiment 7 (solid line with cross). Experiment 8 (solid line) and for 
TMI observations(solid line with dots). 

Fig.8 Time series (3-h intervals) of a) maximum wind (m/s) at the lowest model level 
(about 50 m) and b) minimum sea-level pressure (hPa). 
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